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Abstract 

The decay of A- hypernuclei without n emission, known as Non Mesonic Weak Decay (NMWD), 
gives an effective tool to investigate AS=1 four-baryon interactions, ft was theoretically sug- 
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gested that the two-nucleon induced mechanism could play a substantial role in reproducing 
the observed NMWD decay rates and nucleon spectra, but at present no direct evidence of 
such a mechanism has been obtained. The FINUDA experiment, exploiting the possibility to 
detect both charged and neutral particles coming from the hypernucleus decay, has allowed us 
to deduce the relative weight of the two nucleon induced decay rate to the total NMWD rate. 
The value of T2N /TNMWn=0.24±Q.03 s tat-o < ^2 S s ys s nas been deduced, with an error reduced by 
a factor more than two compared with the previous assessment. 

Key words: A-hypernuclei, non mesonic weak decay, two nucleon induced decay 
PACS: 21.80.+a, 25.80.Pw 



1. Introduction 

The Non Mesonic Weak Decay (NMWD) of Hypernuclei has stimulated a strong in- 
terest since the beginning of Hypcrnuclear Physics [1]. Indeed it has been realized that 
the two channels of NMWD, 

Z A A^ A - 2 (Z-l) + n+p (1) 

and 

f A -> A ~ 2 Z + n + n (2) 

where A and Z indicate, respectively, the mass and atomic numbers of the decaying 
systems, are due to the occurrence of the two weak reactions: 

Ap -> np (T p ) (3) 

and 

An -» nn (T„) (4) 

inside nuclei. These two decay processes are usually referred to as "one-proton induced 
NMWD" and, respectively, "one-neutron induced NMWD" of A-Hypernuclei. Reactions 
(3) and (4) constitute a unique class of four-baryon, strangeness non-conserving weak 
interactions and the determination of their rates is of considerable interest. These rates 
cannot be determined experimentally by the direct reactions (3) and (4), due to the lack 
of suitable beams of Hypcrons. The only way is to study them through their occurrence 
in Hypernuclei (reactions (1) and (2)). 

Unfortunately, also this last approach is not experimentally easy, due to the low rate of 
production of A-Hypernuclei in ground or excited states and to the detection efficiency 
of neutrons, which is small despite the large acceptance and good energy resolution 
apparatuses. For these reasons the experimental progress in this field has been limited 
for many years, in contrast to the vivid development of many theoretical approaches. 
Recent reviews can be found in [2] and [3]. Rcf. [4] first pointed out the possibility that 
a considerable amount of the strength for NMWD, up to about 20% of the total decay 
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width, could be accounted for by the interaction of a A with a pair of correlated nuclcons 
in a nucleus, such as: 



referred to in the following as 2N-induced NMWD. The total non mesonic rate is given 
by rNMWD=ruv + F 2 n, with Tiat= r„ + T p and T 2 n= r„ p + T pp + r„„. 
In addition to reaction (3), (4), (5), (6) and (7), for light hypernuclei one can also study 
the so-called rare two-body decay. Recently FINUDA has been able to measure the decay 
yields and the branching ratios of the two-body decay channels ^He — > d+d, \Ke — > p + t 
and ^Hc -+d + t [5]. 

The suggestion of the existence of the 2N-induced NMWD was followed by detailed cal- 
culations thus stimulating a big experimental effort on the subject. A first summary of 
the theoretical issues can be found in [2], whereas [6] ([7]) reports the latest developments 
in the calculation of the NMWD widths (nucleon spectra) . First experimental determi- 
nations of r2./v were done by indirect methods based on the fit of the experimental 
inclusive proton spectra using IntraNuclear Cascade (INC) calculations including 2N- 
induced NMWD [8] and reported for the T 2 n I^nmwd ratio a value as large as 40%. 
Coincidence measurements of neutrons and protons, with an energy threshold of 30 McV, 
following ^C NMWD were recently analyzed taking into account the angular correlations 
between the detected nucleons [9]. The coincidence spectra were analized using a new 
version of the INC code, with a strength varied to fit the 12 C(p,p') total inelastic cross 
section, in order to account for the effect of the Final State Interaction (FSI) on the ex- 
perimental spectra. From this analysis the experimental value r2Ar/rATMWD=0.29±0.13 
was reported for 12 C [9]. 

A different approach to extract the strength for the 2N-induced NMWD was followed 
by the FINUDA Collaboration [10]. Proton energy spectra of |He, ^Li, ^Be, ^B, x | C, 
^C, ^N and were measured with good resolution (Ap/p=2% FWHM for protons 
of 80 MeV) and with a kinetic energy detection threshold of 15 MeV. All the measured 
spectra showed a similar behaviour, i.e., a bump at about 80 MeV, roughly at the energy 
expected from reaction (1). The bump is quite well defined in the high energy portion, 
whereas at low energies it is blurred in a continuum generated by FSI and superimposed 
to the 2N-induccd NMWD contribution. Under very simple hypotheses, the contributions 
from FSI and 2N-induced NMWD were disentangled, providing: r 2 jv/r p =0.43±0.25 and 
T 2N /Y NMWD =Q.2A±Q.m 

In this Letter we present an improvement of this method, made possible by the detection 
of both protons and neutrons, by which we determine the value of T 2 n /^nmwd with 
an error reduced by a factor of more than two with respect to the previous FINUDA 
assessment [10]. 

2. The experimental and analysis method 

The data were collected by the FINUDA experiment, installed at one of the two in- 
teraction regions of the DA$NE (e + ,e~) ^-factory of Laboratori Nazionali di Frascati 
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(INFN— Italy) and correspond to an integrated luminosity of about 1.2 fb^ 1 . A detailed 
description of the FINUDA experiment can be found in [11,12]. 

We do not report here experimental and analysis details already described in [10] and 
references therein. Some information concerning the performances of the detectors which 
are relevant for the discussion of this new analysis are here recalled together with the 
features of the neutron detection. 

FINUDA was a magnetic spectrometer, immersed in a uniform solenoidal magnetic field 
of 1 T and optimized for the detection of charged particles, with an angular coverage of 
- 2tt sr [11]. The outer FINUDA detector, called TOFONE [10,11] was a barrel of 72 
plastic (CH) n scintillator slabs (255 cm long and 10 cm thick), used essentially for trigger 
and charged particle P.Id. (by Time Of Flight). It was also used to detect neutrons (and 
photons) with an efficiency of about 10% for neutrons in the kinetic energy range 15-150 
McV. 

The inner FINUDA detector was an hodoscope of 12 scintillator thin slabs (TOFINO) 
arranged around the beam pipe at the (e + , e~ ) interaction point; it was used for trigger 
purposes and to identify the charged kaons discriminating them from minimum ionizing 
(Bhabha or beam background) particles. The particle Time-Of-Flight was measured by 
the (TOFINO-TOFONE) system. 

Neutrons and photons are identified looking for events in which TOFONE elements were 
not connected to the curved trajectories belonging to charged particles. 

The /3 of the neutral particle is evaluated by means of j3—(bof/tof)/c, where tof is the 
measured time between the hit, on the TOFINO, of the K~ from the <j> decay just before 
stopping in the target (stopping time<200 ps), and the hit on the TOFONE and &o/is 
the corresponding base of flight. 

The analysis of the /3 values of the neutral candidates allows us to discriminate neutrons 
from 7's emitted, for istance, in the tt° decay. Fig.l shows the distribution of 1//3 for 
selected neutral candidates. The peak centered at 1//3=1 is due to 7's and it is followed 
by the contributions due to neutrons and to 7's from ir° decay, which may be delayed 
if they are produced in the decay of kaons or hyperons. A contamination from charged 
particles back-emitted, following 7 interaction in the iron yoke of the FINUDA magnet 
and from neutron scattering in the FINUDA apparatus, is also present. 
Applying a cut on the 1//3 value (1//3>1.47 as indicated by the line in Fig.l), it is possible 
to eliminate the 7 peak from the neutron candidates detection. As it will be discussed 
in the following the most effective way to identify neutrons from NMWD is to tag them 
when they are accompaincd by a coincidence proton. In the analysis we selected all the 
events in which a neutron and a proton are emitted in coincidence with a 7r _ having a 
momentum corresponding to the A-hypernucleus formation in its ground state or in a 
low lying excited one, decaying to the ground state by electromagnetic emission. 

For the calibration of the neutron energy scale we used the monochromatic neutrons 
(185 MeV/c) produced in the decay at rest of the S + coming from the K~ top + p —> 
S + + 7r~ reaction. The experimental peak is centered at 187.6 MeV/c with a a=9A 
MeV/c. The precision on the determination of the impact point of the neutrons on the 
scintillator slab is tr=6 cm; by combining this value with the timing resolution of tr=780 
ps we finally obtain for the overall energy resolution on the neutron ~13% at 10 MeV 
and ~20% at 100 McV. A detailed description of the TOFONE performances can be 
found in [13]. 

We analyzed neutron spectra coming from NMWD of Hypernuclei by using the same 
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Fig. 1. Distribution of for neutral candidates; the line indicates the cut applied in the data analysis. 

procedures adopted for the proton spectra described in [14], i.e., by requiring the co- 
incidence of a 7r~ from the K~ interaction vertex with a momentum compatible with 
the formation of the Hypernuclcus with a bound A, but the result was unsuccessful. 
The neutron spectra were affected by a huge background, due to the contaminations 
described above, that we could not reduce neither by applying suitable cuts on the ex- 
perimental spectra nor with the help of simulations. Wc then considered the neutron 
spectra obtained by requiring not only the presence of a 7r~, but also of a proton in 
quasi b.t.b. (back-to-back) correlation with the neutron (cos(9(np)<-0.8). The number 
of such triple coincidence events was quite low, typically of the order of twenty for each 
nuclear target, and we could not infer from their distribution reliable conclusions. Wc 
added the events from all hypernuclear species, and we compared neutron and proton 
spectra obtained after the acceptance correction and the subtraction of the background 
due to the K{np) — > Y,~p absorption, followed by the inflight £~ — > ir~n decay [14]. 

Fig. 2 (first row) shows the result; the two spectra, 2a) for protons and 2b) for neu- 
trons, arc quite similar, as expected. The proton coincidence with a quasi b.t.b. neutron 
enhances the number of events due to the process (3) with respect to the 2N-induccd 
channels. On the other hand the quasi b.t.b requirement reduces the number of events 
in the low energy region, which are due to the channel (3) and (4) followed by a re- 
scattering inside the nucleus (FSI); in addition the b.t.b. correlation reduces the number 
of events from the 2N- induced decay (5), which are expected to exhibit a typical 3- 
body phase space angular and energy distribution. A confirmation of such expectation 
was found by requiring a tighter angular correlation (cos#(np)<— 0.9) and the result is 
reported in Fig. 2c) for protons and 2d) for neutrons: the low energy region for both neu- 
tron and proton spectra is depleted and this confirms the validity of the analysis method. 

Taking into account the efficacy of the method described above, we have analyzed the 
triple (7T _ , n, p) coincidences considering the total number of events corresponding to 
each hypernucleus. We showed in [10] that the single proton spectra from NMWD of 
^He and the p-shcll Hypcrnuclei ^Li, ^Be, "B, ^C, ^C, and x \ O had a similar 
behaviour, i.e, a bump around 80 McV, due to reaction (3) without FSI, and a rise in 
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Fig. 2. Upper row: kinetic energy spectra of proton (a) and neutron (b) with a quasi b.t.b. correlation 
(cos9(np)<-0.8) detected in coincidence with a ir~ with a momentum compatible with the formation 
of a A hypcrnuclcus; lower row: kinetic energy spectra of proton (c) and neutron (d) with an angular 
correlation of (cosf?(np)< — 0.9) detected in coincidence with a ir~ with a momentum compatible with 
the formation of a A hypcrnuclcus. 



the low energy region, increasing with A, due to FSI and to channel (5). A fit to each 
spectrum beyond 80 MeV by using a Gaussian function with free mean values and sigma 
allowed us to disentangle the contributions due to process (3), the effect of FSI and 
process (5) and to determine: r 2 Ar/r p =0.43±0.25 and r 2 7v/rAr A /ivD=0.24±0.10 [10]. 
In the analysis of the events due to triple (ir~ , n, p) coincidences first of all we fixed for 
each hypernucleus a proton kinetic energy limit E p placed at 20 MeV below the mean 
value fi of the Gaussian found in [10]. In order to enhance the contribution of the 2N- 
induced NMWD we also chose an upper limit for the proton-neutron angular correlation 
of cos6>(np)=-0.8. 

We classified then the triple coincidence (it~ , n, p) events into four groups, according to 
the following criteria: 

(i) events with proton kinetic energy E p larger than the limit and cos6>(np)< —0.8. 
These events should stem mainly from the process (3) without FSI of the proton, 
even if a negligible contribution from the reaction (5) is expected. Note that we 
neglect the process (6), since T pp ~ 0.15r np [15]; 
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(ii) events with E p larger than the limit and cos#(np)> —0.8. These events should 
correspond mainly to the process (3) or (4) followed by a FSI in addition to a small 
contribution from (5). 
(hi) events with E p lower than the limit and cos#(np)< —0.8. These events should 

correspond mainly to the process (3) with FSI of the proton, 
(iv) events with E p lower than the limit and cos#(np)> —0.8. These events should 
correspond mainly to the process Anp —> nnp with a small contribution of FSI. 
The cuts applied in the selection (iv) produce an underestimation of the population of 
this group which can be estimated to be 20%. This estimation is based on the assumption 
of a >4-body phase space distribution of cos#(np); for the understimation due to the cut 
on E p the calculation of [16] has been used. 

The neutron acceptance for these events was evaluated by taking into account the ap- 
paratus geometry, the efficiency of the FINUDA pattern recognition algorithm and the 
quality cuts applied to the real data. The emission of a neutron from a K~ stopped in 
the target was simulated, assigning to the neutron, emitted in coincidence with a pro- 
ton and a negative pion, a fiat momentum distribution from 100 MeV/c to 700 McV/c. 
The proton momentum was also simulated with a flat distribution in the same range of 
the neutron and the 7r~ momentum was simulated in the range (270-290) MeV/c corre- 
sponding to the formation of the A hypernuclei. The acceptance function was evaluated 
target by target in different apparatus sectors and it was found to be flat within 20% in 
the energy range (10-100) MeV. 

3. Experimental results and discussion 

We considered the dependence on the mass number A of the ratio R between the 
events of the selection (iv) and the number of protons with E p >n. This last number 
was evaluated by integrating the number of events of the proton spectra above the mean 
value /i of the Gaussian fit of [10] for each hypcrnuclcus. We have: 

_ N n (E p <(fi- 20MeV),cos6(np) > -0.8) _ 
= N P {E P > n) ~ 

0.8N{Anp -> nnp) + N^ SIlN + N^ sl2N 
0.5N(Ap np) + Np SIlN 

where N n {E p < (a — 20 M eV), cos 9 (np) >-0.8) corresponds to the number of neutrons 
fulfilling the condition (iv) and that we estimate to be 80% of the total number of neutron 
coming from the process (5), N p (E p > n) is the number of protons with an energy higher 
than the mean value of the fit reported in [10]; N£ SXlN , N^ SIlN , N^ ShN arc the number 
of protons or neutrons which suffer the final state interaction effect following the one- 
(1N) or two-nucleon (2N) induced NMWD. In Fig. 3 the experimental values of this ratio 
for each hypernucleus are plotted as a function of A. 

Due to FSI both N n (E p < (n - 20MeV), cos6(np) > -0.8) and N p (E p > u) in (8) 
are expected to be proportional to A; we performed a fit with a function R(A) = (a' + 
b'A)/(c + dA) and we found d/c<10~ 3 ; so a simple linear fit R(A) = (a + bA) can be 
applied. The result of such a fit (a = 0.63 ± 0.25, b = -0.017 ± 0.021 and x 2 /ndf = 
0.298/6) is reported in Fig. 3. 

We can make the assumption that r np /r p is constant in the mass number range of the 
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Fig. 3. The ratio N n (E p < Q-t - 20MeV), cos0(np) > -0.8)/N p (E p > fi) as a function of the mass 
number A. 

analyzed hypernuclci [6]. Moreover, one has: 



thus eq.(8) provides: 



N(Anp -» nnp) _ T^p 
N{Ap ->np) ~ r p ' U 

where b is taken from the fit and represents the small FSI contribution. Eq.(lO) can be 
solved for T np /T p for each hypernucleus, giving: 



IV = [R(A) - bA] 
T p 1.6 



(11) 



The values obtained for T np /T p for the examined hypernuclei are all compatible with 
each other within errors. The final result can thus be given by the weighted average over 
all the available nuclear species: 

r 



= 0.39±0.07 stat ™;- . (12) 

1 jn 



!!£_nqnj-nn7 +0.04 3 
p 

The quoted statistical error is less than half the one obtained in the proton spectra 
analysis of [10]. This is because, unlike [10], using the described selection criteria, the 
ratios R and T np /T p are linearly related. The experimental inaccuracies due to the back- 
ground subtraction and the efficiency correction are included in the statistical error. The 
systematic errors instead take into account the effect of slight variations in the selection 
criteria applied on cos#(np) and on the limit of the proton energy E p ; the variation due 
to the use of the weighted average instead of the normal one is also included. In addi- 
tion, the possible uncertainty on the understimation of population of group (iv) has been 
included. 

As pointed out in [10], the 2N-induced NMWD can be assumed to be dominated by 
the Anp —¥ nnp channel and the recent microscopical calculation of [15] delivers T np : 
r pp : r„ n = 0.83 : 0.12 : 0.04. Ref. [15] presents an improvement of the first calculation 
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of the 2N- induced decay rates performed within a nuclear matter framework [4] , where a 
phenomcnological description of the two-particle two-hole polarization propagator was 
adopted; the effects of Pauli exchange terms in the two-nucleon stimulated NMWD are 
also taken into account in [15]. Considering that the total error on r np /r p quoted in 
eq.(12) is comparable with the contribution of (T pp + T nn ), in the present analysis we 
cannot assume T2N ~ r np as in [10]. We can determine I^tv/Tnmwd by adopting the 
same method employed in [10]; using the experimental value of T n /T p reported in [8] for 
A He and ^C, our determination of T np /T p and taking into account that T np /T nmwd — 
0-S3(T2n/^nmwd) [15] we obtain: 

f ~~ (r 7r \ 4-1 I cr — Ttm ~ u U.03 s t at _ 02 . (13) 

iNMWD U n/ 1 p) + J- + (1 2N/Lp) v 

This value supports the latest theoretical predictions [6] (r 2 /rNMWD=0-26), the recent 
experimental results of [9] (0.29 ±0.13) and the previous FINUDA result [10], but bears 
a smaller error. 



4. Conclusions 



We performed a model-independent analysis of the 2N-induced NMWD process Anp — > 
nnp by analyzing the (7r _ , p, n) triple coincidence events from K~ stopped in thin 
nuclear targets with the FINUDA spectrometer at DA$NE. The measurement of the 
7r~ momentum allowed the selection of events coming from the decay of p-shell A- 
hypernuclei. By applying appropriate cuts on the energies and angles of the protons and 
the neutrons we could identify events due mainly to Anp — > nnp NMWD and determine 
r„ p /r p =0.39±0.07 sta ^o;Q3^ s s . We also extracted r2N/rNAiWD=0-24±0m s tattl' o ll v y l re- 
ducing the total error by a factor of more than two with respect to the value published in 
[10]. The obtained results arc in agreement, within the errors, with previous experimental 
determinations, model dependent or not, and with theoretical calculations. We hope that 
these limited statistics results will be improved by the experiments planned at J-PARC. 

We dedicate this paper to the memory of Dr. Ambrogio Pantalco passed away prema- 
turely on June 23rd, 2010. Dr. Pantaleo participated very actively since the beginning to 
the FINUDA experiment, taking particular care of the neutron detection, on which the 
present analysis is based. 
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